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ABSTRACT 

Pseudolites can extend the availability of GPS-type 
positioning systems to a wide range of applications not 
possible with satellite-only GPS, including indoor and 
deep-space applications. Conventional GPS pseudolite 
arrays require that the devices be pre-cali bra ted through a 
survey of their locations, typically to sub-centimeter 
accuracy. This can sometimes be a difficult task, 
especially in remote or hazardous environments. By 
using the GPS signals that the pseudolites broadcast, 
however, it is possible to have the array self-survey its 
own relative locations, creating a Self-Calibrating 
Pseudolite Array (SCPA). 


In order to provide the bi-directional ranging signals 
between devices necessary for array self-calibration, 
pseudolite transceivers must be used. The basic 
principles behind the use of transceivers to create an 
SCPA were first presented in f 1 ]. This paper begins with 
a brief review of the transceiver architecture and the 
fundamental direct-ranging algorithm presented in that 
paper. This is followed by a description of a prototype 
self-differencing transceiver system that has been 
constructed, and a presentation of experimental code- and 
carrier-phase ranging data obtained using that system. A 
second algorithm is then described which uses these 
fundamental range measurements between transceiver 
pairs to self-calibrate a larger stationary array and to 
provide positioning information for a vehicle moving 
within that array. Simulation results validating the 
accuracy and effective convergence of this algorithm are 
also presented. 

INTRODUCTION 

The robotic exploration of other planets poses many 
challenges for navigation systems. Familiar terrestrial 
position and orientation references may be limited. Man- 
made landmarks such as roads do not exist, and neither do 
radio navigation aids such as GPS, LORAN, and 
VOR/DME. Moreover, the long communications time 
delay and limited bandwidth limit human intervention and 
dictate a requirement for a great deal of system autonomy. 

Mars surface exploration suffers from many of these 
navigational challenges. The lack of a global magnetic 
field makes heading determination difficult without the 
use of star trackers. Poor traction prevents the use of 
odometers for long-range positioning, and rugged terrain 
and long-duration missions limit the use of intertia! 
instrument packages. 

These difficulties are accompanied by a growing demand 
for accurate and reliable navigation information. Sample- 
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return missions require positioning accuracy at the sub- 
meter level, in order to allow robots to locate and return to 
valuable scientific sites. As mission plans transition to 
more distributed systems with multiple cooperating 
robots, even more accurate relative positioning 
information is needed to enable successful completion of 
mission goals. This requirement also exists for potential 
astronaut/robot teams in the future. 

GPS-type systems can greatly benefit future space 
exploration missions by providing centimeter-level, drift 
free position and degree-level attitude information to 
robotic explorers. GPS pseudolites allow this technology 
to be used in areas where satellite coverage is not yet 
available, or to supplement the satellites with additional 
ranging sources. The basic feasibility of this approach 
has been previously demonstrated. For example, the use 
of GPS pseudolite augmentation for relative positioning 
of formation-flying spacecraft is described in f 21~f7] . 
Ground-based augmentation for situations with poor 
satellite coverage is described in f 8]-[ 1 0]. 



Figure I : Mars SCPA 


Although the motivation for this research is Mars 
exploration, SCPAs can be used in other terrestrial and 
space applications where the surveying of the pseudolites 
by other methods is difficult, impractical, or otherwise 
impossible. The algorithms and experimental results 
presented in this paper are general, and are applicable to 
most such applications. 


Pseudolite-only systems have also been previously 
demonstrated (e.g. for indoor positioning systems fill- 
113]). A similar pseudolite system, distributed on the 
Martian surface, would allow precise positioning of 
multiple robots over extended duration missions. A major 
difficulty with such systems, however, is in surveying the 
locations of the pseudolites: these must be known with an 
accuracy comparable to the overall navigational accuracy 
desired from the array. This is a significant challenge for 
autonomous rovers in a hostile environment. Zimmerman 
demonstrated that it is possible to survey in a pseudolite 
system using their broadcast signals by moving a receiver 
within the array to precisely known locations, thereby 
inverting the pseudolite navigation problem [13]. 
Although this method is infeasible for planetary 
exploration, the ability to have the array self-calibrate 
using the broadcast GPS signal is still highly desirable. 


TRANSCEIVER ARCHITECTURE 

Transceivers can be constructed using a variety of 
architectures. Many of the possibilities are summarized 
in [14]. The architecture chosen for this project is a self- 
differencing transceiver, with separate transmitter and 
receiver components, as shown in Figure 2. This 
architecture was chosen for ease of construction and 
component integration. The output of the pseudolite is 
split, with one line going to a passive broadcast antenna 
and the other going to one front end on the dual front-end 
receiver. This allows the receiver to monitor the output 
signal, effectively measuring the relative clock bias 
between the pseudolite and its own receiver. The other 
receiver RF front end is connected to an antenna that 
listens for signals from other transceivers. 


Such self-calibration is possible by collocating GPS 
receivers with each pseudolite transmitter, thereby 
creating a GPS transceiver. An array of transceivers is 
able to listen to the signals broadcast by each member of 
the array and determine the overall configuration, creating 
a Self-Calibrating Pseudolite Array (SCPA). Positioning 
to a couple of meters accuracy is possible using code- 
phase measurements, while centimeter-level accuracy is 
possible using carrier-phase measurements once the 
integer ambiguity is resolved. Figure 1 shows a 
conceptual drawing of such a system providing 
navigational capability to a rover on the Martian surface. 
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INTER-TRANSCEIVER RANGING 


The simplest navigation solution using sell-differencing 
transceivers directly determines the range between the 
antennas on a pair of devices, using only the signals 
broadcast and received by the devices themselves. This is 
also a very useful solution, because the ranges between 
any number of devices can all be solved for individually 
and with no changes to the fundamental equations. Figure 
3 shows such a pair of devices. A description of the 
terminology used appears below. Note that the time 

biases 7 ; and 7 J include common-mode effects such as 
transmit line biases up to the splitting point of the signal. 
The measurement 0 t J can represent either a carrier- or 

code-phase measurement. In the case of carrier-phase 
measurements, the integers can simply be included with 

the line biases b / , and either solved for or removed by 
using a-priori position information. 

bj Line bias from PL j to Rec i 
Ry Range between device antennas 

(p- Rec i’s measurement of PL j 
7 ( . Clock bias of Rec i 
7^ Clock bias of PL j 



Figure 3: Inter-Transceiver Ranging 


For this analysis it is assumed that the airborne RF 
pathlength between each of the two transmit/receive 
antenna pairs is identical and equals the range between 
the devices. The antenna geometry presented in Figure 3 
satisfies this constraint, although other geometries may be 
used with slight changes to the equations. It is also 
assumed that the measurements are latched 
simultaneously. The errors associated with violating this 
assumption are described in [I] and [13], and can 
generally be made negligible over short (< 1000 km) 
baselines by using intelligent receiver latching strategies. 

The raw measurements taken by each receiver of the 
signals from the two pseudolites are given in Equation 1. 
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Eliminate any receiver clock biases or common mode 
effects by taking internal single (self) differences between 
the signals received by a given receiver, as shown in 
Equation 2. 

a<*. = « >/ -<#; = (*/ (r 1 

Wi **i -«>; =4; -*;)+ (r J 


Combining the measurements from both receivers and 
inverting Equation 2, one can determine both the range 
between the antenna pairs and the relative clock bias of 
the pseudolites. 
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Note that unlike most conventional differencing schemes, 
the line biases do not cancel. This is because the signals 
do not travel along common lines. A similar effect occurs 
in multiple-antenna attitude systems. These biases must 
therefore be removed, either by a hardware calibration or 
by solving for them explicitly in the navigation 
algorithms. 

Noise characteristics of the range measurement between a 
pair of transceivers are the same as for a conventional 
single-difference measurement, i.e. twice that of the raw 
receiver code- or carrier-phase measurement noise. 

For very long baseline measurements (tens to thousands 
of km) two additional factors become important. First, as 
was mentioned earlier, drift in the receiver and pseudofite 
clocks make it difficult to latch the measurements 
simultaneously. This effect can be alleviated somewhat 
by using higher-grade oscillators. A second effect is that 
unlike conventional double-differencing schemes, 
ionospheric and troposheric delays present in terrestrial 
applications do not cancel. Rather, they average and add 
directly to the measured range. 
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EXPERIMENTAL SYSTEM 

This section describes the prototype system used to test 
and validate the inter-transceiver ranging technique 
presented in the previous section. The GPS components 
of the system are two self-differencing transceivers. The 
pseudolite portion of each transceiver is an IntegriNautics 
IN200C signal generator, shown in Figure 4, broadcasting 
differing pseudorandom codes. These generators were 
chosen because of their many desirable features: 
programmable power levels, several different pulsing 
schemes, a programmable data message at variable data 
rates, and the ability to broadcast at frequencies offset 
from LI. Figure 5 shows the receiver used; a Mitel Orion 
receiver that has been modified to include two separate 
RF front ends [15]. These receivers are fully 
programmable, allowing the modification of the tracking 
loops to accommodate the pseudolite data structure. Each 
receiver is equipped with an RS-232 serial link for data 
collection. 



Figure 4: Pseudolite Signal Generator 



Figure 5: GPS Receiver 


1 he antenna configuration used for these transceivers is 
shown in Figure 6. The upper black antenna is a passive 


patch antenna used for broadcasting the pseudolite signal, 
and the lower white antenna is an active patch antenna 
that receives signals from other transceivers. The vertical 
circular groundpianes help reduce multipath interference 
from behind the antennas and also add directionality to 
the antenna patterns. The horizontal groundplane 
between the antennas reduces the transmission between 
the transmit and receive antennas of the same device, 
helping to alleviate the near-far problem. 



Figure 6: Antenna Configuration 


Code- and carrier-phase data were collected 
simultaneously from both transceivers in the pair at a rate 
of 10 Hz, and was then post-processed to determine the 
range between the devices. Each data packet was time 
stamped by the receiver based on the current code -epoch, 
allowing the data to be synchronized. Missing data points 
were filled in by interpolation and cycle slips were 
removed. Real-time data processing was also performed, 
and showed similar performance, but did not insure 
synchronicity of the data or correct for missing data 
points. 

EXPERIMENTAL RESULTS 

The experimental data in this section come from a scries 
of tests of the prototype system described above in an 
indoor laboratory setting. The baselines between the 
transceiver pair were relatively short (0.5 to 3.0 m), 
giving a difference in signal strength between near and far 
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operations of about 15.6 dB. This is smaller than the 
difference encountered in long-baseline operations (the 
collocated pseudolite is broadcasting at a very high signal 
strength), but allows the use of the Orion receivers 
without any modifications to the tracking loops. Future 
long-baseline (-100 m) tests are planned, following 
receiver modifications to allow tracking of very high SNR 
signals. During these tests the pseudolites were 
generating standard C/A signals; no pulsing was used. 

The indoor test environment is a relatively cluttered room 
about 8 m wide and 9 m long, containing a significant 
amount of equipment. It therefore presents a severe 
multipath environment. No special efforts were taken to 
reduce the effect of this multipath, which is believed to be 
worse than the system would see in outdoor operations. 
Because of the excellent results obtained from the indoor 
system, major multipath problems are not anticipated 
outdoors. 

Although the data collection and processing system 
checks for missing data points and corrects for cycle slips, 
these have shown themselves to be infrequent events. 
The receivers skip data points every 1-2 hours. Carrier- 
phase cycle slips occur with similar frequency for slow 
(<20 cm/s) relative motion between the transceivers, 
although they become more frequent at faster relative 
motions and at receiver SNRs below 10-12 dB. Code- 
phase cycle slips occur more frequently (several times per 
hour), but are correctable in a similar manner. Moreover, 
their large magnitude makes them easily detectable. 

Figure 7 presents carrier-phase ranging data from the 
indoor test rig. The antennas started at a range of 3.0 m, 
and were then moved inward to a closest range of 0.5 m 
in 0.5 m increments. Truth measurements were provided 
by a metric scale on the floor, and are accurate to 
approximately 1-2 cm at antenna height. It was assumed 
that the initial starting separation (and thus the integer 
ambiguity and line biases) was known. Figure 8 shows 
the corresponding SNR values of the received signals in 
one of the two receivers, both from its own collocated 
pseudolite and from the other transceiver. The near-far 
change in SNR is clearly evident. 

The ranging data show excellent tracking of the carrier- 
phase throughout the course of the test and through a wide 
range of SNR values. The mean positioning error for all 
of the stationary placements together was 1.29 cm, about 
the level of the technical error in the truth system. This 
includes the 4.7 cm error for the Final placement, which is 
most likely due to operator error (the system placement is 
discretized at 5 cm). The other errors are due to a 
combination of the uncertainty in the true antenna 
locations and other error sources such as multipath. The 
ranging measurement noise level is very low, with a 
standard deviation of 2.3 mm. 


Carrier-Phase Ranging 




Figure 9 shows code-phase ranging from similar test. The 
devices were initially separated by 3.0 m, were moved to 
1.0 m, and finally back to 3.0 m. The upper plot shows 
the raw code-phase ranging data, while the lower shows 
carrier-smoothed code using a complementary filter 


R + ; Carrier 


s + a 


s + a 


(4) 


which is discretized using the bilinear transformation. 
The break frequency a = 0.001 rad j s forces very close 
tracking of the carrier-phase data, but allows a steady- 
state offset due to the steady-state code-phase range 
estimate. The break frequency can be chosen with this 
low value because there is no intervening ionosphere to 
generate code and carrier divergence. For this plot the 
initial range estimate is the average value of the code- 
phase range over the first two minutes of data, corrected 
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Code-Phase Spectrum 


for the line-biases estimated from another collection of 
data using the same hardware setup. 


Code-Phase Ranging (m) 
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Figure 9: Code-Phase Data 

Although the code-phase data track the motion of the 
devices well, the roughly half-meter offset in the data 
shows the difficulty in accurately determining the code- 
phase line biases. Typical runs may have steady-state 
offsets from each other as large as 1-2 m. In addition, any 
given data set will see very low frequency variations of 
this magnitude as well. Figure 10 shows a smoothed 
empirical transfer-function estimate (ETFE)[I6] of the 
frequency spectrum of 30 minutes of raw (unsmoothed) 
code-phase data, clearly exhibiting non-white noise 
characteristics with a more significant low-frequency 
component. This low-frequency variation is difficult to 
remove on a real-time basis, and the effect can be 
considered similar to a low-amplitude S/A-type error 
source. Very long averaging windows would help 
alleviate this problem in situations where there is no 
relative motion between the transceiver pairs for long 
periods of time. 


SELF-CALIBRATION ALGORITHMS 

SCPA self-calibration can be done using either code- or 
carrier-phase measurements. Code-phase measurements 
give limited (meter-level) accuracy, but have the 
advantage of not having an integer ambiguity. This 
means that once line biases have been calibrated out, a 
single set of static measurements serves to completely 
self-survey the array. For greater accuracy carrier-phase 
measurements are used, although this adds the 
aforementioned ambiguity. Motion of one of the 
transceivers can serve to resolve this ambiguity, and is the 
approach taken for the initialization algorithm presented 
here. 



The algorithm presented in this paper differs from 
previous formulations because it uses the computed range 
between transceiver pairs - as determined using the intcr- 
transceiver ranging technique described above - as the 
fundamental measurement type. It therefore presents a 2- 


step calculation process, first finding the 



ranges between the N transceivers in the array, and then 
combining these range measurements to determine the 
array positions. Depending on the size and 
dimensionality of the array some of these range 
measurements may be redundant; these can either be 
incorporated using a least-squares type method, or some 
can simply be eliminated and ignored. 


This method of using inter-transceiver ranges as an 
intermediate step provides several advantages over other 
algorithms. Recordkeeping is simplified, because all 
possible permutations between transmit and receive 
antennas do not need to be considered. Range is a 
fundamental physical quantity and thus makes system 
evaluation and diagnostics simpler, and also makes it 
easier to determine which measurements are redundant 
and can be eliminated. This greatly helps simplification 
when the array includes a large number of devices. 
Finally, the 2-step process reduces the number of 
unknowns that must be determined at once, offering a 
potential increase in algorithm speed. 


Other differencing schemes involving more than two 
transceivers as the fundamental unit are less intuitive and 
more difficult to implement, but can offer potential 
advantages. One such advantage is the ability to use a 
mobile receiver on its own, instead of a full transceiver, to 
resolve integer ambiguities. This reduces equipment 
requirements at the expense of modularity and simplicity. 
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Because of the near-field nature of the SCPA, the problem 
of finding the device locations is inherently non-linear. 
The algorithm presented here uses standard linearized 
iterative least-squares to move from a poor initial estimate 
of array and vehicle position to a final accurate estimate. 
This initial estimate can be provided by direct code-phase 
ranging between the transceivers or through other sensors 
including intertial instruments and computer vision 
systems. 

The nomenclature used in this discussion is summarized 
below.. . 

b Vector of range measurement biases 

R n Range measurements at step /? 

£ n Measurement errors at step n 

P n Vehicle position at step n 
V Vector positions of static transceivers 

These values at a given vehicle position are consolidated 
into common vectors and matrices to solve as a batch 
process. 


ih - 


dGjx) 

dx 


a/?, a/?. 


-i 


dv 


dv ap, 

a^, 


o 


0 


a^ 

dp. 


( 7 ) 


Because this only approximates the non-linear relation, an 
iterative approach is used to descend the gradient. At 
each iteration step k, the solution is given by computing 

the left pseudoinverse H 1 . 


$r h =n' w k = (u T k nJn k dv, w 


where 


= v m 




(9) 
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The carrier-phase measurements H 7 collected over the 
course of the trajectory are related to the state vector X 
through the non-linear relation 


is the difference between the actual and expected 
measurement values. This provides the estimate of the 
state which minimizes the difference between the actual 
and the expected measurement values in a least squares 
sense. This effectively minimizes the following cost 

function J k . 

J, = J (‘W, - H t &, ) ( 10 ) 


V = G(x) + E 


(5) 


The state update equations are then 


To solve this relation, it is convenient to linearize the 
system using small variations from the assumed state 

values x k . 

=H t Sc k +(5E* (6) 


where 


** + . =** -& k 
=G(jf w ) 


(ID 


Note that in order to compute the pseudoinverse // 1 one 
must constrain the solution. This can either be done by 
defining the coordinate system based on the transceiver 
locations and eliminating the corresponding columns of 

H k , or by adding constraint equations by a technique 

such as Lagrange multipliers. Because of its simplicity 
the former method is adopted here, although it does skew 
the estimation error slightly towards one side of the array. 
One transceiver is defined to be at the origin and another 
is defined to lie along the x-axis. This provides the 

constraints necessary to make II k full rank and 
invcrtable. 
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SELF-CALIBRATION SIMULATION 


In order to verify that this algorithm accurately converges 
to the actual system states using a realistic initial estimate 
of transceiver locations and rover path, the following 
simulation was performed. The hypothetical test situation 
is a stationary, planar triangular array of three transceivers 
with 100 m separation. A fourth transceiver is mounted 
on a vehicle that is free to move about within the array. 
This is the minimum number of devices required to 
resolve all the observables for this case; there are no 
redundant measurements. It is assumed that carrier-phase 
measurements are being taken, and that the vehicle 
motion will be used to resolve both the integers and the 
uncalibrated line biases. 


To examine the effect of geometry change on the 
determination of the integers and line biases, four 
different trajectories of the vehicle through the array were 
examined: a straight line through the middle of the array, 
a lawnmower pattern in the middle of the array, a 
complete circuit around the outside of the array, and a 
figure-eight type pattern which crosses through the array 
between loops around the static transceivers. These four 
patterns are shown in Figure 11. In each case, 
measurements are taken at roughly 10 m intervals. 


T rajectory A 
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Figure 11: Sample Vehicle Trajectories 

A singular value decomposition (SVD) of the 

pseudoinverse of the linearized geometry matrix fl L 
was used to evaluate the relative effectiveness of each of 
these trajectories. The error in the final state estimates is 
proportional to the measurement errors divided by the 

singular values O r 


Trajectory C 




< 7 . 


( 12 ) 


Values of the minimum singular value O min for each 
trajectory are presented in Table 1. Singular values close 
to or greater than one will give cenlimcter-Icvcl position 
estimates, and are considered very good for evaluation 
purposes. Trajectory A performed quite poorly. 
Trajectory B was slightly better, but still did not provide 
adequate observability of all the states. Trajectories C 
and D both performed extremely well, with excellent 
observability of all states. This is because the looping 
around the outside of all three static transceivers. Tests 
show that eliminating the loop around even one 
transceiver can reduce the observability by a factor of 10 
or greater. This is similar to the effect driving a vehicle 
around a pseudolite to resolve the integer ambiguities 
associated with CDGPS, as presented in f 17]. 


Table 1: Minimum Singular Values for Trajectories 


Trajectory 

A 

B 

C 

D 

O 

mm 

0.0038 

0.0172 

0.6087 

0.5977 


Because of its excellent resolution and because it presents 
an interesting path, Trajectory D was used for the 
algorithmic simulation. The unconstrained degrees-of 
freedom of the static transceivers were given a random 
initial estimate. The vehicle was assumed to start at a 
location and with a heading close to its estimated starling 
point, and then proceeded to attempt to follow the 
assigned trajectory. 

Tracking errors included both random variations in the 
distances traveled and the steering angle and small linear 
drift terms in both as well. These errors are similar to 
those exhibited in odometers when wheel slippage occurs 
and when steering angles are not known precisely, and so 
have some level of physical plausibility. They were also 
chosen because they present a more difficult challenge to 
the algorithm, which is able to cope with spatially random 
initial position and trajectory estimates of up to half the 
size of the array. The magnitude of these error 
parameters, which are listed in Table 2, were chosen to 
give a realistic trajectory tracking accuracy which could 
be achievable with other non-GPS navigational sensors. 
In addition, unknown biases (integers and line biases) 
with a standard deviation of 1000 m were added to each 
range measurement. Measurement errors w'ere selected to 
be typical of CDGPS systems, with a 1 cm standard 
deviation. 


Table 2: Trajectory Tracking Prior Parameters 



Bias 

Linear Drill 

Random 

Distance 

10 m 

0. 1 m/m 

0.1 m/m 

Angle 

5° 

0.0257m 

0.57m 
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Figure 12 shows both the initial guess of the trajectory 
and transceiver locations and the final estimates. Both the 
trajectory and the transceiver locations converged well, 
giving an RMS error of 0.0074 m over the entire array. 
Figure 13 shows the RMS error in the estimated states as 
the iteration process continues. In this case, the errors fall 
below one centimeter RMS after 12 iterations. Bias 
estimate errors fall to the same order of magnitude as the 
other estimate errors in one iteration. Multiple simulation 
runs with different random trajectories show that this 
algorithm is reasonably robust, and performs well (80- 
90% successful convergence) with estimation parameter 
errors of twice the magnitude presented here, even though 
the actual trajectory at this point bears little resemblance 
to the initial path estimate. 

Continuing algorithmic work is focused on using code- 
phase ranging, potentially blended with other sensor data, 
to provide more accurate initial estimates for the carrier- 
phase iteration process. Other algorithms are also being 
examined to attempt to reduce the sensitivity to the initial 
state estimate. 



CONCLUSIONS 

This research demonstrates the potential benefits of direct 
ranging between transceivers for positioning applications, 
and demonstrates the accuracy of a self-differencing 
transceiver architecture for this purpose. The data 
presented show inter-transceiver ranging accuracy 
comparable to that of standard CDGPS. Of particular 
interest is the utilization of inter-transceiver ranging as the 
fundamental measurement unit for the surveying of 
SCPAs. The simulation shows successful array 
calibration from a reasonable initial state estimate, using a 
mobile transceiver to resolve integers. Future field tests 
will further validate the SCPA concept with full hardware 
demonstrations. 
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